What the Gribov copy tells on the confinement and 
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We performed lattice Landau gauge QCD simulation on (3 — 6.0, 16*, 24*, 32* and (3 = 6.4, 32*, 48* 
and 56* by adopting the gauge fixing that minimizes the norm of the gauge field, and measured the 
running coupling by using the gluon propagator and the ghost propagator. In view of ambiguity 
in the vertex renormalization factor Zi in the lattice, we adjust the normalization of the running 
coupling by the perturbative QCD results near the highest momentum point. It has a maximum 
Qs(g) ~ 2.1(3) at around q = 0.5 GeV and decreases as q approaches 0, and the Kugo-Ojima 
parameter reached -0.83(2). The infrared exponent of the ghost propagator at 0.4GeV region is 
ac = 0.20 but there is an exceptional Gribov copy with ac = 0.27. The features of the exceptional 
Gribov copy are investigated by measuring four one-dimensional Fourier transform(l-d FT) of the 
giuon propagator transverse to each lattice axis. We observe, in general, correlation between absolute 
value of the Kugo-Ojima parameter and the degree of reflection positivity violation in the 1-d FT of 
the gluon propagator. ©The 1-d FT of the exceptional Gribov copy has an axis whose sample- wise 
gluon propagator manifestly violates reflection positivity, and the average of the Cartan subalgebra 
components of the Kugo-Ojima parameter along this axis is consistent to -1. The running coupling 
of the ensemble average shows a suppression at momentum, but when the ghost propagator of the 
exceptional Gribov copy is adopted, the suppression disappears and the data implies presence of the 
infrared flxed point Qs(0) ~ 2.5(5) and k = 0.5 suggested by the Dyson-Schwinger approach in the 
multiplicative renormalizable scheme. Comparison with the SU(2) QCD and Nf = 2 unquenched 
SU(3) QCD are also made. 

PACS numbers: 12.38.Gc, 11. 15. Ha, ll.15.Tk 



I. INTRODUCTION 

The lattice Landau gauge QCD simulation suffers from 
Gribov copy problem and its effect on the confinement 
was discussed by several authors Q, |^ |M Q . As a method 
for obtaining the unique gauge, we adopted the funda- 
mental modular gauge (FMG) i.e. a configuration with 
the minimum norm of the gaiige field and studied the 
Gribov copy problem in SU(2) 5]. We compared the ab- 
solute minimum configuration obtained by the Landau 
gauge fixing via the parallel tempering method and the 
1st copy which is obtained by our straightforward Lan- 
dau gauge fixing. We observed that the FMG configu- 
rations and the 1st copy which is in the Gribov region 
but not necessarily in the FM region have the following 
differences: 1) The absolute value of the Kugo-Ojima pa- 
rameter c 1^ , which gives the sufficient condition of 
the confinement, of the FMG is smaller than that of the 
1st copy. 2) The singularity of the ghost propagator of 
the FMG is less than that of the 1st copy. 3) The gluon 
propagator of the two copies are almost the same within 
statistical errors. 4) The horizon function deviation pa- 
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rameter h of the FMG is not closer to 0, i.e. the value 
expected in the continuum limit, than that of the 1st 
copy. 

The proximity of the FMG configuration and the 
boundary of the Gribov region in SU(2) in 8*, 12* and 
16* lattices with [i = 0,0.8,1.6 and 2.7 was studied in 
8J. The tendency that the smallest eigenvalue of the 
Faddeev-Popov matrix of the FMG and that of the 1st 
copy come closer as (3 and lattice size become larger was 
observed, although as remarked in Q the physical vol- 
ume oi (3 — 2.7, 16* lattice is small and not close to the 
continuum limit. Qualitatve features of the profile of the 
Morse function 

m = \Y. I d^^{[A(^^T{x)Y (1) 

where g = e*^ "^, was sketched as a function with respect 
to the magnitude of the infinitesimal gauge transforma- 
tion parameter e and a parameter r which is defined by 
2nd, 3rd and 4th derivative of £[g\ with respect to e at 
the origin. The simulation suggests that as the (3 and lat- 
tice size become large, the parameter r decreases. The 
meaning of the parameter r is such that larger r than 
the critical value implies an existence of a smaller local 
minimum than that of the origin. 

The difference of the 1st copy and the FMG in the 
(3 = 2.2, 16* lattice W\ indicates that the FMG does 
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not overlap with the boundary of the Gribov region in 
that simulation. In the Langevin formulation of QCD, 
Zwanziger conjectures that the path integral over the FM 
region will become equivalent to that over the Gribov re- 
gion in the continuum 3J. This conjecture is consistent 
with the view that the boundary of the FMG and that 
of the Gribov region overlaps and the probability distri- 
bution is accumulated in this overlapped region. On the 
lattice, when /3 and the lattice size is not large enough, 
distribution of Gribov copies i.e. statistical weight of the 
copies is crucial for extracting sample averages. 

In the previous paper Q, we measured the QCD run- 
ning coupling and the Kugo-Ojima parameter in P — 
6.0, 16", 24", 32" and (3 = 6.4,32"* and 48"*. The running 
coupling was found maximum of about 1.1 at around 
q = 0.5 GeV, and behaved either approaching constant 
or even decreasing as q approaches zero, and the Kugo- 
Ojima parameter was getting larger but staying around 
—0.8 in contrast to the expected value —1 in the con- 
tinuum theory. Thus it is necessary to perform a larger 
lattice simulation and to study the dependence of the 
Gribov copy. We encountered a rather exceptional Gri- 
bov copy in P — 6.4, 56'* which is close to the Gribov 
boundary and we consider it worthwhile to investigate 
that sample in some details. We analyze those data by 
comparing with continuum theory like Dyson-Schwinger 
equation(DSE). 

There are extensive reviews on DSE for the Yang Mills 
theory 0, IHIll El- The solution of DSE depends 
on ansatz of momentum truncation and what kind of 
loop diag rams are included. Two decades ago Mandel- 
stam [13 projected the DSE for the gluon propagator by 
'Pfiiyil) = — qiiqvl'f' and without including ghosts, 
assumed the gluon wavefunction renormalization factor 
in the form 

^(9') = -T + ^(9') h = const. (2) 
T 

Later Brown and Pennington jl4| argued that in order 
to decouple divergent tadpole contribution, it is more 
appropriate to project the gluon propagator by TZ^,^(q) = 
5 — Aq^qi, / q^ . A careful study of inclusion of ghost loop 
in this DSE was performed by [ll|, and they showed the 
infrared QCD running coupling in Landau gauge could 
be finite. 

The divergent QCD running coupling caused difhculty 
in the model building of dynamical chiral symmetry 
breaking In order to get reasonable values of the 

quark condensates, infrared finite QCD running coupling 
was favored. Recent DSE approach with multiplicative 
renormalizable(MR) truncation with infrared finite QCD 
running coupling llM ITflj l suggests that the confinement 
and the chiral symmetry breaking can be explained by 
the unique running coupling. We thus compare the run- 
ning coupling obtained from our lattice simulation and 
that used in the DSE and study the dependence on the 
Gribov copy. 

We produced SU(3) gauge configurations by using the 



heat-bath method, performed gauge fixing and analyzed 
lattice Landau gauge configurations of /? = 6.4, 56". The 
(3 = 6.4,48" and 56" lattices allow measuring the ghost 
propagator in the momentum range [0.48,14.6] GeV, and 
[0.41,14.6] GeV, respectively. In the present work, the 
gauge field is defined from the link variables as log U 
type: 

TJ — 4t — — A 

The fundamental modular gauge (FMG)|2| of lattice 
size L is specified by the global minimum along the gauge 
orbits, i.e., 

Al = {U\Fu{l) = MingFuig)}, Al C ^l, where 
is called the Gribov region (local minima) and 
= {U\ - dD{U) > , dA{U) = 0}. 
Here Fu{g) is defined as 

In the gauge transformation 

gA«.,^ = gle'^^-'^gx+f,, (3) 

where g = e*^'^, the value e is chosen depending on the 
maximum norm |9yl|cr as follows. 

. When \dA\ > \dA\,r: = p^^^- iv' - 0.05) 

• When \dA\ < \dA\^r- e = {-df,D^{A))-^T]dA (77 = 
1 - 1.6) 

In the second case, calculation of (— 9^I?^(v4))~* is 
performed by Newton's method where the linear equation 
is solved up to third order of the gauge field, and then 
the Poisson equation is solved by the multigrid method 
The accuracy of the gauge fixed configura- 
tion characterized by dA(U) = is 10^" in the maximum 
norm squared which turned out to be about 10^^^ in the 
L2 norm squared of the gauge field in contrast to about 
10-"2 in 48". 

In the calculation of the ghost propagator, i.e. inverse 
Faddeev-Popov (FP) operator, we adopt the conjugate 
gradient (CG) method, whose accuracy of the solution in 
the q < 0.8GeV region turned out to be less than 5% in 
the maximum norm j^, Q ■ 

In , we analyzed these data using a method inspired 
by the principle of minimal sensitivity (PMS) and/or the 
effective charge method the contour- improved 
perturbation method (2^ and the DSE approach [lol[T5| . 
We perform the same analysis to the 56" data. 

The infrared behavior of the running coupling is tightly 
related to the mechanism of the dynamical chiral symme- 
try breaking 17, 18, 27]. The lattice data are compared 
with the theory of dynamical chiral symmetry breaking 
based on the DSE. 

In order to study properties of Zi and the infrared 
features, we extend the 16" SU(2) lattice Landau gauge 
simulation and compare data oi (3 = 2.2, 2.3 and 2.375. 
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In sec. II we show some details of the gauge fixing pro- 
cedure and show sample dependence of the gluon prop- 
agator, Kugo-Ojima parameter and QCD running cou- 
pling. In sec. Ill a brief summary of the DSE as well as 
the recent exact renormalization group approach(ERGE) 
are presented. We compare lattice data with results 
of the theoretical analysis of DSE. The SU(2) lattice 
Landau gauge simulation data are summarized in sec. 
IV. In order to check qualitative differences between the 
quenched and unquenched Landau gauge simulation, we 
performed an exploratory analysis of the configuration 
produced by the JLQCD^J- The results are shown in 
sec. V. Summary and issues on dynamical chiral symme- 
try breaking is discussed in sec. VI. 

II. GRIBOV COPY AND THE 56* LATTICE 
DATA 

The magnitude of |i9^|cr in the gauge transformation 
is chosen to be 2.2(copy A) or 2(copy B). In most cases, 
gauge fixed configurations are almost the some, but in 
some cases, different |9A|cr produce significantly differ- 
ent copies. 

In order to see the difference of the gluon field of 
the Gribov copies, we measured the 4 components of 1- 
dimcnsional Fourier transform (1-d FT) of the sample- 
wise gluon propagator as follows. We consider the gluon 
propagator 

DaMi) = tr{A,,{q)A,{qy) 

^ iS,^^^)DA{q% (4) 

where A^{q) ~ —j= e^*^^A^(x). In the data analysis 

ofim there are some possible choice of q. Here we choose 
q transverse to [i. Since there are 3 possible choices of 
z/ 7^ /z, we make an average of the three combinations 

= (5) 

When the axis v is chosen as t axis, and an average 
over is taken, it is equivalent to the specific Schwinger 
function 

1 

S{t, 6) = — Y, DAiqo, 0)e2-««*/^ (6) 
qo=0 

where L is the lattice size. 

When the Schwinger function becomes negative, the 
reflection positivity becomes violated, which means that 
the gluon is not a physical particle. Violation of pos- 
itivity is considered as a sufficient condition of the 
confinement [Toll27Ll2^. 



The four 1-d FT of the copy Ia and those of the copy 
Ib are shown in Fig. ^and in Fig. 13 respectively. The 
solid line, dotted line, dashed line and the dash-dotted 
line corresponds to propagator transverse to xi,X2,X3 
and X4 axis in the Euclidean space, respectively. 

lOOi 
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FIG. 1: The 1-d FT of the gluon propagator along the 4 axes. 
/3 = 6.4,56* in the log U definition, sample Ia 
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FIG. 2: The 1-d FT of the gluon propagator transverse to the 
4 axes. f3 = 6.4,56* in the log U dfinition. sample Ib 

We observe that the gluon propagators of copies Ia 
and Ib have a specific axis along which the propagator 
manifestly violates reflection positivity. Here, manifestly 
means that it remains negative in a wide range in the 
intermediate not only in the large distance in the coordi- 
nate space. Propagators transverse to other axes in the 
copy Ib are shifted from those of Ia and the propaga- 
tor almost parallel to that manifestly violating reflection 
positivity remains finite in the copy I a, but it becomes 
almost in the large distance in the copy Ib- The L2 
norm squared of copy Ib is smaller than that of Ia, 
and hence Ib is closer to the FMR but is not necessar- 
ily closer to the boundary of the Gribov region. Rather 
small shifts of the gluon propagators among copies make 
a significant difference in the exponent of the ghost prop- 
agator and the Kugo-Ojima parameter is surprising. 
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The ghost propagator is defined by the expectation 
value of the inverse Faddeev-Popov(FP) operator M 

D'i{x,y) = {{\^x\{M[U])-^\XS)), (7) 

via the Fourier transform 

Doi^) = (8) 

The Kugo-Ojima parameter is defined by the two point 
function of the covariant derivative of the ghost and the 
commutator of the antighost and gauge field 

= ^ E ^"^'^'"'^ (^"^^M 4^ [Au , A"]) ) . (9) 

We performed the same analyses as sample / for a 
sample which has the second largest Kugo-Ojima param- 
eter (samples Ha and IIb)- The sample dependences of 
the L2 norm of the gauge field, Kugo-Ojima parameter 
c — — m(0), trace divided by the dimension e/d, horizon 
function deviation parameter h 0, l20l | and the infrared 
exponent of the ghost propagator at 0.4GeV region ac, 
are summarized in Table |I| Errors in c is due to the de- 
viation of the tensor sturucture from — ^^^|^) i.e. c 
depends on the choice of n as in the exceptional copy. 

We parametrize infrared power dependence of DA{q^) 
as - ((7a)-2(i+"^) andthatofDG('7^) as - {qa)~^'-^+°"^\ 
Errors in aa are estimated from the standard deviation in 
the plot of log Dg as a function of log q and we find they 
are (-0.1,-1-0.45). An analysis of DSE '3^ suggests that 
the exponent at 0.4GeV is about half of the asymptotic 
value K and thus ac corresponds to about half of k. From 
our standard deviation of aa, we expect k in the range 
of [0.1,0.7]. 

TABLE I: The Gribov copy dependence of the Kugo-Ojima 
parameter c, trace divided by the dimension e/d, horizon con- 
dition deviation parameter h and the exponent ac- 



sample average is 0.22, but that of the Ia copy is 0.27. 
The I A copy has a larger L2 norm of the gauge field but 
smaller h and larger c. We find that not all samples have 
the axis that manifestly violates reflection positivity and 
that the direction of the axis is sample dependent. 

A. Kugo-Ojima parameter 

Our sample average of c = — m(0), e/d, h, the exponent 
of the ghost dressing function ac, the exponents of the 
gluon dressing function an near q = 0.4Gey, and 
near q — 1.97GeV are summarized in Table ITTl 

TABLE IL The Kugo-Ojima parameter c, trace divided by 
the dimension e/d, horizon function deviation h in the log U 
definitions. The exponent of the ghost dressing function near 
zero momentum ac, the exponent of the gluon dressing func- 
tion near zero momentum q_d, near q = 1.97GeV a'^ in logi7 
type. f3 = 6.0 and 6.4. 





6.0 


6.4 


L 


16 


24 


32 


32 


48 


56 


c 


0.628(94) 0.774(76) 0.777(46) 


0.700(42) 0.793(61) 0.827(27) 


e/d 


0.943(1) 


0.944(1) 


0.944(1) 


0.953(1) 


0.954(1) 


0.954(1) 


h 


-0.32(9) 


-0.17(8) 


-0.16(5) 


-0.25(4) 


-0.16(6) 


-0.12(3) 


aa 


0.175 


0.175 


0.174 


0.174 


0.193 


0.223 


ao 




-0.310 


-0.375 




-0.273 


-0.323 


ol'd 


0.38 


0.314 


0.302 


0.31 


0.288 


0.275 



The color off-diagonal, space diagonal part of the 
Kugo-Ojima parameter c was 0.0001(162) and consistent 
to 0. The magnitude of the Kugo-Ojima parameter c and 
exponent of the ghost propagator ac are tightly corre- 
lated and they are also correlated with the violation of 
the reflection positivity in the gluon propagator. In the 
I A copy, reflection positivity is violated along axis and 
the average of 33 and 88 color components of c along this 
axis is 0.97(6), consistent with 1. 





Ia 


Ib 


IIa 


IIb 


average 


\\Af 


0.09081 


0.09079 


0.090698 


0.090695 


0.09072(7) 


c 


0.851(77) 


0.837(58) 


0.835(53) 


0.829(56) 


0.827(15) 


e/d 


0.9535(1) 


0.9535(1) 


0.9535(1) 


0.9535(1) 


0.954(1) 


h 


-0.102(77) 


-0.117(58) 


-0.118(53) 


-0.125(56) 


-0.127(15) 


aa 


0.272 


0.241 


0.223 


0.221 


0.223 



We observed that in most samples the dependence of 
the copy on |9A|cr is weak as in the case of sample //, 
and that the large difference of Ia and Ib copies is ex- 
ceptional. The Table m also shows that ac c and h are 
correlated. In the average of 15 samples of 56^ lattice 
data, we incorporate copy A but not B. The ac of the 



B. Gluon propagator 

The gluon propagator in momentum space was mea- 
sured by using cylindrical cut method l29l |. i.e., choosing 
momenta close to the diagonal direction. In Fig. O we 
show the gluon dressing function of /? = 6.4, 56** lattice 
data together with 48''^ lattice data. The gluon propa- 
gators of 24^,32^ and 48^ as a function of the physical 
momentum agree with each other within errors and they 

can be fitted by the MOM scheme in two loop pertur- 
bation theory 0, ■ 

DA{q') = ^^^''^^'r°-°'^'^ = (10) 
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The overall normalization in this fitting turned out to 
be problematic since the 56^ data are suppressed than 
the 48^ data. We remove the lattice artefact by rescaling 
the data of the dressing function to that of the fit in the 

MOM scheme ZA{9.5GeV) = 1.3107(9) 





FIG. 4: The ghost propagator as the function of the momen- 
tum g(GeV). 13 = 6.0, 24* (star), 32* (unfilled diamond),^ = 
6.4, 48* (triangle) and 56* (filled diamond) in the log?7 defini- 
tion. The fitted line is that of the MOM scheme. 



FIG. 3: The gluon dressing function as the function of the mo- 
mentum g(GeV). /3 = 6.4, 48* (triangles) and 56* (diamonds) 
in the log [7 definition, extrapolated to V = oo. The solid 

line is that of the MOM scheme. All data are scaled at 
IJ, = 9.5GeV. 



C. Ghost propagator 

The ghost dressing function is defined by the ghost 
propagator as G"^''(g2) = q^Da'^^q^)- In Fig. H /3 = 6.4, 
48"*, and 56"* and fi = 6.0 24* and 32* lattice data of the 
ghost propagator are compared with that of the MOM 
scheme |^l30j|. 



Zg{q^,y)\v=o. 



02142 



(11) 



We observe that the agreement is good for q > 0.5 



GeV. The MOM scheme is singular at Aj;jg ~ 0.35 GeV 
but the singularity should be shifted to momentum 
by the non-perturbative effects. The ghost propagator 
was first measured in [sJl but the scaling property was 
not observed and the lowest momentum point was incor- 
rectly suppressed. It may worth while to remark that the 
rescaling is unnecesary in the ghost propagator of differ- 
ent lattice sizes, but the scale depends on the definition 
of the gauge field. The propagator of log U definition is 
about 14% suppressed from that of the U— linear defini- 
tion. 



D. QCD running coupling 

We measured the running coupling from the product of 
the gluon dressing function and the ghost dressing func- 
tion squared p^l44|. In terms of exponents ao and ac 



the running coupling near 0.4GeV is parametrized as 



as{q') 



ZA{q^)G{q^ 



(qa) 



-2(aD+2QG) 



(12) 



The lattice size dependences of the exponents ao and 
ac are summarized in Table Hll 

The vertex renormalization factor Zi is 1 in the per- 
turbation theory, but on the lattice it is not necessarilly 
the case. By comparing data of various 0, finiteness of Zi 
was confirmed in the case of SU(2) In the present 

analysis, we fix Zi by normalizating the running cou- 
pling by that of the perturbative QCD near the highest 
momentum point. In the lattice simulation of the three 
gluon coupling ^40], the nonperturbative effect is found 
to be significant even at 10 GeV region, and a fit of the 
lattice data by the three loop perturbative term plus c/q^ 
correction was proposed. We normalize the running cou- 
pling to that of Orsay group at the point of 14.4 GeV, i.e. 
0.154(1). This correction revises the previous results of 
48* lattice data0 by a factor of 1.97, and the maximum 
of the running coupling becomes 2.0(3). 

In Fig[Sl we present the rescaled running coupling of 
48* lattice and that of the 56* lattice and the fit of Orsay 
group above 2GeV and the result of the MR truncation 
scheme of BlochJ^ _19,], where in addition to the sunset 
diagram, the squint diagram was included. The running 
coupling in this DSE is parametrized as 



oisiq 



1 



Co 



47r 
(Jo logt 



t - 1 



)f) (13) 



where t — q^/Ag^jj. The infrared fixed point ag is ex- 
pressed as an analytic function of k, and ^3 claims that 
when two-loop squint diagrams are included, possible so- 
lutions exist only for k in the range [0.17,0.53]. Conjec- 
tures from DSE [iTL IT^ predicts k ~ 0.5, which implies 
ao - 2.5. 



6 



Except the value of the lowest momentum point, our 
data is consistent with the prediction ao = 2.5. Thus, we 
adopt this value for ao and search the parameter cq by the 
fit to the second and the third lowest momentum points 
of the running coupling. We find parameter cq = 30, 
instead of cq = 15 in the DSE [l3|. 

Phenomenologically fitted Aqcd from a{Mz) is about 
710 MeV, but the value depends on the number of quark 
flavors and in the quenched approximation the choice is 
not appropriate. We choose as |18| . 



Aqcd = 330 MeV. 




2 4 6 8 



10 



12 14 



FIG. 5: The running coupling as(q) as a function of momen- 
tum g(GeV) of the /3 = 6.4, 56** lattice and 48* lattice. The 
DSE approach with ao = 2.5 (long dashed line) and the the 
Orsay group(dotted line) are also plotted. 

When the ghost propagator of the exceptional copy is 
adopted, suppression of the running coupling at mo- 
mentum disappears. The DSE results, Orsay fit and the 
lattice data of the running coupling in which the ghost 
dressing function is taken from the average as a func- 
tion of logarithm of momentum logj^g '7(G'eV) are shown 
in Figl^l In order to show the dependence on the Gri- 
bov copy, the data in which the ghost dressing function 
is replaced by that of the Ia copy is also shown in the 
same figure. The ensemble of gluon propagator was not 
changed in this replacement, since the sample- wise dif- 
ference of the gluon propagator is insignificant. 

The contour improved perturbation method with A = 
g70/6/3o^__ jppp Qpjej- jg consistent with 

our data at q > lOGeV region, (dotted hue) but in the 
infrared region it underestimates the lattice data. The 
dotted line is qualitatively the same as the results of hy- 
pothetical T lepton decay |35l |. 



III. COMPARISON WITH DSE AND ERGE 

In the DSE approaches, infrared power behavior and 
specific relation between the exponent of the ghost prop- 
agator and the gluon propagator is assumed. In the 
ERGE, flow equation in terms of the effective aver- 
age action Fa where A is the infrared cut-off scale is 



4r 

3.5- 

3- 
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5- 
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0.5- 
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FIG. 6: The running coupling as{q) as a function of the log- 
arithm of momentum logjp g(GeV) of the f3 = 6.4, 56'* lat- 
tice using the ghost propagator of the Ia copy (stars) and 
that of the average (diamonds). The DSE approach with 
ao ~ 2.5(long dashed line), the fit of the Orsay group pertur- 
bative +c/q'^ (short dashed line) and the contour improved 
perturbation method ( dotted line) are also shown. 



considered |32l l33l l34j. In a recent work four point ver- 
tices in addition to the two point vertices are incorpo- 
rated and the running coupling was calculated via 



47r/2(g2;A^0)/G'(<z2;A^0) 



(14) 



where fz{q^', A) and fciq^'-, A) are gluon and ghost prop- 
agator function, respectively. They are related to the 
gluon and ghost propagator as 

DA,t,i.{q^) = - %qu/q^)^^-n-^, — — (15) 



and 



Daiq^) = - 



1 



g2/G(g2;A^0) 



(16) 



The infrared exponent k obtained in this analysis turned 
out to be K ~ 0.146 in contrast to the DSE approach 
which suggested k ^ 0.5. The infrared fixed point 
ao ^ 4.70 was predicted [l^ which is about factor 2 larger 
than our lattice simulation. There is a prediction k — 
0.59535 • • • and the infrared fixed point ao = 2.9717- • • 
both in DSE and ERGEJ2? ]. 

The prediction ao = 2.6 and k — 0.5 of is con- 
sistent with our lattice data. Here we summarize his 
approach and compare our lattice results. 

The quark prop agator in Euclidean momentum state 
is expressed as [l^l27| 



1 



Z{q') 



(17) 



and M{q'^) = B{q^) / A{q^) is proportional to the quark 
condensate at large g^: 



M(g2) 



A.TTas{q^) ^as{q^) . 



3^2 



'q;s(m^ 



'(#(m')) (18) 
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where dm = 12/(33 — 2Nf). Here the number of flavor 
Nf = in the quenched approximation. 
The quark field is renormalized as 

Z{q^,^^^) = Z2{^l^A^)ZR{q^^i^) (19) 

where Zji is the renormaUzed quark dressing function, 
Z2 is the quark field renormalization constant and at the 
renormahzation point. We define Zn^x) = Zr(x, p!^) and 
ZniiM^) = 1 and = Min"^). In the DSEfia, /x is 
chosen to be Aqcd = 330MeV. 

The renormaUzed quark dressing function Zji{q) and 
the quark mass function M{q) can be calculated by a cou- 
pled equation once the running coupling as{q^) is given 
[l8l| . The quark mass function at the origin M(0) is a 
function of the parameter Cq and our fitted value cg = 30 
yields 

M(0) ~ 1.27Aqcd = 0.419GeV. (20) 

This value is consistent with the result of quark propaga- 
tor in quenched lattice Landau gauge simulation 33 ex- 
trapolated to momentum. The quark condensate ('0V') 
is estimated as — (0.70Aqci5)'^ which is compatible with 
the recent analysis of quenched lattice QCD [38l |. 

When K is larger than 0.5 as predicted by 0? ], 
the gluon propagator should vanish in the infrared. The 
present lattice data are not compatible with this predic- 
tion. 

IV. SU(2) 16" LATTICE DATA 

In finiteness of the vertex renormalization fac- 

tor Zi was prooved by hnear rising of the Z a{p.^)G{ii'^)'^ 
{pL = 3GeV) as a function of — \og{a{(3)'^ a) where a(/3) 
is the lattice spacing corresponding to the /3 and a = 
[AAQM eVY' is the string tension. In order to check this 
behavior and to see infrared features of the SU(2) lattice 
Landau gauge, we performed Monte Carlo simulation of 
SU(2) lattice Landau gauge using the J7— linear defini- 
tion of the gauge field. We choose (3 — 2.2, 2.3 and 2.375 
and accumulate 200 samples for each (3. 

We confirmed increasing of Z a{h'^)G{im^Y (m = 3GeV) 
from (3 = 2.3 to 2.375 with the slope 7 consistent with 
13/22. The data of /3 = 2.2 was oflF the fitted line, but 
we expect this is due to the closeness of the ^ = iGeV 
point to the maximal momentum point 3.7GeV. 

The gluon dressing function and the ghost dressing 
function as a function of the momentum g(GeV) of the 
(3 = 2.2, 2.3 and 2.375 are shown in FigEland in FiglHl 
respectively. In the gluon dressing function, cylindrical 
cut is applied and the error bars are obtained by the jack- 
nife method. Error bars of the ghost dressing function 
are the standard deviation. 

The running coupling as{q) as a function of the log- 
arithm of the momentum log2o[9(GeV)] oi 13 — 2.2,2.3 
and 2.375 are plotted in FigO We normalize the run- 
ning coupling near the highest momemtum point by that 
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FIG. 7: The SU(2) gluon dressing function as a function of 
momentum q(GeV) of the P = 2.2(triangles), 2.3(diamonds) 
and 2.375(stars), 16* lattice (200 samples). 
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FIG. 8: The SU(2) ghost dressing function as a function of 
momentum g(GeV) of the (3 = 2.2(triangles), 2.3(diamonds) 
and 2.375(stars), 16^ lattice (200 samples). 



of the two- loop perturbation results. This correction re- 
vises the previous result of the running coupling of SU(2) 
J] by about factor 1.54, but there remains difference from 
Tuebingen and Sao Carlos '4^ by about factor 2. 

The ghost propagator and the gluon propagator of |43 | 
were rescaled by the tadpole renormalization factor up. 
There are qualitative agreements in ghost propagator of 
[4^ 1 and ours, but in the gluon propagator there are dis- 
crepancies in the momentum dependence in the infrared 
region. In [45l | it is remarked that in simulations of rela- 
tively small lattice with a lattice axis chosen to be twice 
as those of the other three lattices, the gluon propagator 
of a few lowest momentum points do not match smoothly 
to those of higher momenta. Our gluon propagator be- 
low 1 GeV is more suppressed than those of Tuebingen 
data of 16"^ x 32J14J, and the discrepancy could be due 
to this finite size effect. Tuebingen group adopts adjoint 
links and thus their tadpole renormalization makes di- 
rect comparison of the gluon propagators obscure. In 
the running coupling, however, the tadpole renormaliza- 
tion factor u p for the ghost and for the gluon cancel j 
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FIG. 9: The SU(2) running coupling as{q) as a function 
of the logarithm of the momentum log]^Q[g(GeV)] of the 
f3 = 2.2(triangles), 2.3(diamonds) and 2.375(stars), 16'' lat- 
tice (200 samples). The result of two- loop perturbation theory 
(dotted line ) are also plotted. 



and the difference in the running coupling as{q) in the 
infrared is due to the difference in the shape of the gluon 
propagator. 



V. UNQUENCHED SU(3) 20^ x 48 LATTICE 
DATA 

We observed that there are samples whose 1-d FT of 
the gluon propagator transverse to a lattice axis mani- 
festly violates reflection positivity. The direction of the 
reflection positivity violating axis appears randomly. Re- 
cently, Aubin and Ogilvie 46] pointed out that the origin 
of the reflection positivity violation lies in the quenched 
character of the gauge transformation g. They demon- 
strared in a Higgs model type SU(2) 20'' lattice simula- 
tion, occurlence of reflection positivity violation analo- 
gous to that in the quenched lattice simulation of oq me- 
son propagator, by considering the gauge transformation 
of Giocai X G global- In order to see qualitative difference 
between quenched and unquanched simulation, and to 
investigate finite size effects in lattices whose one axis is 
taken longer than the others, we studied infrared features 
of the unquenched SU(3) 20^ x 48 lattice configuration 
of the JLQCD,47,], where improved Wilson action with 
Sheikholeslami and Wohlert parameter csw = 2.02 and 
the number of sea quark flavours Nf = 2 are adopted. 
We choose SHMC algorithm conflgurations of hopping 
parameter Kgea — 0.1340 and 0.1355. 

We performed the Landau gauge flxing on 9 samples for 
each Ksea using the log[/ definition for the gauge field 
and measured the gluon propagator, ghost propagator, 
QCD running coupling and the Kugo-Ojima parameter. 
In addition to the correlation of gauge fields around the 
diagonal [gi, 92, qs, 94] = [q, q, q, (48(?/20)] where (48(7/20) 
is an integer close to this quotient, we measured the cor- 



relation transverse to the coordinate 

X—X,t 

= (4. - '-^)DA{q'). (21) 

where k and / run over 1,2 and 3 ^ i, and the same ex- 
pression for the time axis X4. The four Id-FT of these 
sample- wise gluon propagators turned out to be quite dif- 
ferent from those of the quenched simulation. As shown 
in Figs^ and 12 in quenched case in general, there is a 
component which remains positive in the whole region 
and different from other three components, but such a 
component is absent in the unquenched case in general. 
Although symmetry violation in each sample does not 
mean symmetry violation in the ensemble, the differ- 
ence suggests that the global symmetry, i.e. rotational 
symmetry is recovered by the coupling of the gluon to 
fermions [i^ . 

When the length of an axis is taken longer than the 
other three axes, the Z{A) symmetry corresponding to 
interchange of the axes is broken to .Z(3) symmetry and 
it is serious in the estimation of the infrared gluon prop- 
agator. In the 20'^ x 48 SU(3) unquenched lattice simu- 
lation, the gluon propagator of the 3 lowest momentum 
points (pi,P2,P3,P4) =(0,0,0,0), (0,0,0,1) and (0,0,0,2) do 
not match smoothly to higher momenta (Fig llO|) . 

The problem due to lack of rotational symmetry in 
the gluon propagator is usually evaded by performing 
the cylindrical cut. In this context, the momentum 
points (1,0,0,0) and its Z{3) partners are farther than 
the (0,0,0,2) to the cylidrical axis and the treatment of 
these points remains a problem. In a preliminary calcula- 
tion of the running coupling using (1,0,0,0) and its Z(,i) 
partners of the unquenched 20'^ x 48 lattice, we found 
Q!s(0.98GeV)~ 3 and the infrared fixed point of ao ~ 4 
is suggested. The gluon propagator obtained by Landau 
gauge fixing unquenched SU(3) configuration in which 
Liischer-Weisz improved action is used shows that the 
rotational symmetry of 20'^ x 64 lattice is recovered[? ]. 
Details of the investigation of the running coupling of 
unquenched SU(3) Landai gauge simulation will be pre- 
sented elsewhere. 



VI. DISCUSSION AND OUTLOOK 

We measured the gluon dressing function and the ghost 
dressing function in lattice Landau gauge QCD and cal- 
culated the running coupling. In view of uncertainty in 
the vertex renormalization factor Zi which is not neces- 
sarily 1 in the lattice simulation, we normalized the run- 
ning coupling by that of perturbative QCD near the high- 
est momentum point of the lattice. We found infrared 
fixed point ao ~ 2.5(5), which is consistent with the MR 
scheme DSE calculation [l^. In the momentum depen- 
dence, there is disagreement with DSE in 2 < g < 10 
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FIG. 10: The gluon propagator DA{q) as a function of the 
momentum g(GeV) of the /3 — 5.2, 20^ x 48 lattice using the 
configuration of Ksea =0.1355. 

GeV region which suggests a correction like c/q'^ term 
in as{q) |4ll |. Although this correction applies only in 
q > 2 GeV region, it could yield attraction between col- 
ored sources. 

We observed that the 1-d FT of gluon propagator of the 
Ia copy has an axis along which the reflection positivity 
is manifestly violated. The average of Cartan subalgebra 
components of Kugo-Ojima parameter along this specific 
axis becomes consistent with c = 1. The 1-d FT of the 
gluon propagator transverse to the diagonal direction in 
the lattice is also performed by using the analytical ex- 
pression of the gluon dressing function in MOM scheme 
for q > I GeV and numerical interpolation for < g < 1 
GeV. In this ensemble average, violation of reflection pos- 
itivity is very weak, although the quantitative feature is 
sensitive to the dressing function near q = 0. 

When the QCD running coupling in the infrared re- 
gion is thought to be divergent, the dynamical chi- 
ral symmetry breaking was thought to be irrelevant to 
confinement [l7|. Our lattice data of running coupling 
is qualitatively similar to that assumed in the model of 
dynamical chiral symmetry breaking. 

In passing, we compare running coupling measured in 
other lattice simulations. Orsay group measured the run- 
ning coupling with use of linear definition and from 
triple gluon vertex. The running coupling turned out 



to behave as oc in the infrared contrary to ours, but 
above 0.8 GeV the data are consistent with ours. They 
analyzed the infrared behavior in the instanton liquid 
model '40'|. Running coupling around 0.2 GeV in instan- 
ton scheme u sing U— linear definition measured by the 
DESY group [i^] is = 4 ~ 5. A comparison of the 
logC/ and [/—linear definitions of 48"' and 56** are pre- 
sented in [20I [4^ . The ghost propagator in [/—linear 
definition is larger than that of log U definition, but the 
dependence does not explain the discrepancies in the run- 
ning coupling from the DESY data, and we suspect prob- 
lems in finite size effects due to asymmetric shape of the 
lattice. 

In the study of instantons, Nahm conjectured that Gri- 
bov copies cannot tell much about confinement ,4]. We 
showed that the dynamical chiral symmetry breaking and 
confinement can be explained by using the same running 
coupling and that the Gribov copy gives information on 
the ambiguity in the parameter that characterizes chiral 
symmetry breaking and confinement. We are currently 
analyzing infrared properties of the unquenched JLQGD 
configurations, i.e. the quark mass dependence of the 
Kugo-Ojima parameter and the running coupling. The 
results will be published in the future. 

The running coupling of the quenched SU(3) simula- 
tion of our Landau gauge fixing suggests that there is a 
peak of as ^ 2.2 at q ^ 0.5GeV, but the running cou- 
pling calculated by the ghost propagator of the excep- 
tional sample is consistent with the result of DSE with 
infrared fixed point as{0) ~ 2.5. Whether the population 
of the exceptional configuration becomes larger when the 
system approaches to the continuum limit will be inves- 
tigated. 
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